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When light travels through strongly scattering media with optical gain, the synergy between diffu¬ 
sive transport and stimulated emission can lead to lasing action. Below the threshold pump power, 
the emission spectrum is smooth and consistent from shot-to-shot. Above the lasing threshold, the 
spectrum of emitted light becomes spiky and shows strong huctuations from shot-to-shot. Recent 
experiments have reported that emitted intensity resembles a power-law distribution {i.e. Levy 
statistics). Recent theories have described the emergence of Levy statistics as an intrinsic property 
of lasing in random media. To separate intrinsic intensity huctuations from the motion of scatterers, 
we compare the statistics of samples with stationary or freely-diffusing scatterers. Consistent with 
previous reports, we observe Levy-like statistics when intensity data are pooled across an ensemble 
of scatterer conhgurations. For fixed scatterers, we find exponential intensity distributions whose 
mean intensities vary widely across wavelengths. Levy-like statistics re-emerges when data are com¬ 
bined across many lasing modes. Additionally, we find strong correlations of lasing peak intensities 
across wavelengths. A simple mean-field statistical model captures the observed one- and two-point 
statistics, where correlations in emission intensity arise from competition among all lasing modes 
for limited gain. 


Mesoscopic transport and Anderson localization have 
been widely studied for electrons, photons and cold 
atoms. Unlike electrons and atoms, photons may be 
multiplied via stimulated emission [T]. The synergy of 
coherent amplification and multiple scattering leads to 
fascinating phenomena such as lasing in random media 
las]. From a basic physics point of view, random lasers 
are complex, open, nonlinear systems that bridge various 
fields such as mesoscopic physics, nonlinear dynamics, 
laser physics, and quantum optics. Random lasers ex¬ 
hibit unique characteristics, such as tunable/low spatial 
coherence mi] and spectral fingerprinting of the random 
structure [ 7 ] , which points to a wide range of applications 
in high-speed full-field imaging [8] , optical tagging lano], 
and cancerous tissue mapping min]. 

Over the past two decades there have been extensive 
experimental and theoretical studies on random lasers 
m- One fundamental difference from conventional lasers 
is strong intrinsic fluctuations in the number of lasing 
modes, lasing frequencies, and emission intensities [14]- 
[26] . Such fluctuations are attributed to the existence 
of numerous modes with similar lasing thresholds and 
strong mode interactions via the gain material. One 
interesting observation is that the distribution of emis¬ 
sion intensities can exhibit features of Levy statistics, 
including a power-law decay. The emergence of Levy 
statistics is predicted by analytic and numerical models 
that account for the statistics of random walks but ig¬ 
nore the interference effect iniiiiisi. However, most 
experiments on statistics of random-laser emission were 
conducted with mobile scatterers, with random configu¬ 
ration changes from one pump pulse (shot) to the next 


m UHlEolEIlES]. In general, the strong fluctuation 
of emission spectra from shot-to-shot can arise from (i) 
amplification of noise from spontaneous emission or fluc¬ 
tuations of the pump, (ii) changes in the positions of mo¬ 
bile scatterers. The separate contributions from different 
sources to the random laser statistics are not known. 

In this Letter, we separate the contributions of scat¬ 
terer motion from intrinsic optical fluctuations by com¬ 
paring the intensity statistics of random laser systems 
with identical optical properties but different scatterer 
motion. The data from samples with stationary scat¬ 
terers suggest that individual lasing modes have an ex¬ 
ponential distribution of energies while Levy statistics 
emerges when data are pooled across spectrally hetero¬ 
geneous modes. Additionally, we find strong intensity 
correlations between lasing modes and provide an empir¬ 
ical model that captures the observed one- and two-point 
statistics. 

We consider emission from colloidal scatterers dis¬ 
persed in dye-doped liquids. Polystyrene scatterers have 
a radius of 125 nm and are suspended in aqueous media 
at a number density of about 0.25 /im“^ (0.2 % volume 
fraction). Using Mie scattering theory, we calculate the 
scattering mean free path of these samples to be about 
410 jam. The scatterers are suspended in an aqueous 
solution of 4.5 mM sulfarhodamine 640 and lauramine 
oxide, a surfactant used to enhance the solubility of the 
dye (Ammonyx LO diluted 4:1 with water). 

Samples are optically pumped with 30 ps pulses of 532 
nm light at a repetition rate of 10 Hz from a frequency- 
doubled pulsed NdiYAG laser (Continuum Leopard). 
The pump light is focused by a 2.5 cm focal length lens 
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to a ^ 10/im spot. The axial position of the cuvette is 
chosen to maximize the contrast between laser emission 
peaks and smooth background, which occurs when the 
excitation pulse is focused a few hundred microns into 
the sample. The emitted light in the backward direction 
is collected by the same lens, passing through a long-pass 
filter (to remove the scattered pump light) and then di¬ 
rected to a grating spectrometer (Acton 300i). The exit 
port of the spectrometer is connected to a CCD array 
detector, which is synchronized with the pump laser to 
record the emission spectra of individual pulses. When 
the pump pulse energy exceeds a threshold, discrete las¬ 
ing peaks appear in the emission spectrum. To ensure 
that the lasing peaks are well-resolved above the ampli¬ 
fied spontaneous emission (ASE) background, but also to 
minimize bleaching of the dye, the incident pump pulse 
energy is set to 400 nJ, which is roughly twice the lasing 
threshold. For each sample we record single-shot emis¬ 
sion spectra over 500 pulses at a fixed pump level. 

Results for scatterers immersed in the liquid gain 
medium are shown in the first row of Fig. Raw 

emission spectra, /(A,t), from 200 consecutive shots are 
shown in Fig. [^. The emission spectrum of each shot is 
composed of about twenty discrete peaks, and the peak 
frequencies change completely from shot-to-shot. Thus, 
each peak appears as an isolated bright spot on the 2D 
plot of emission spectra across shots in Fig. [^a). While 
emission peaks dominate single-shot spectra, the emis¬ 
sion spectrum averaged over 500 shots (/(A))^ is smooth, 
as shown in Fig. Bb)- Each wavelength shows large fluc¬ 
tuations in the intensity of emission, and we see no cor¬ 
relation in the location or height of the peaks from shot- 
to-shot. We quantity the intensity fluctuations at each 
wavelength using the survival function, 5'(/, A), in Fig. 
[^c). The survival function, 5'(/, A), gives the fraction 
of shots with emission intensity at wavelength A above 
the intensity /. Compared to density histograms, this 
method of quantifying fluctuations has better fidelity in 
the tail because there is no binning and every intensity 
value is plotted. Note that exponential and power-law 
probability distributions give exponential and power-law 
survival functions. Here, the survival functions fall off 
with a power-law tail with exponent v ^ —1 before reach¬ 
ing a cut-off at near 10 times the mean intensity. This 
is consistent with previous reports m na uni HU [25] 
of Levy statistics of the emission of random lasers. The 
fourth panel shows the distribution of intensities pooled 
across wavelengths. For the red curve, S{I/ the 

intensities are pooled directly across shots and wave¬ 
lengths, and normalized by the mean over all shots and 
wavelengths (J)^ For the black curve, S{I{X)/ (/(A))^), 
the intensities at each wavelength are first normalized by 
the average intensity of that wavelength (/(A))^ before 
pooling across wavelengths. These two distributions are 
essentially identical, suggesting that the distribution of 
intensities is similar across the emission spectrum. 


In the above experiments, the structure of the scat¬ 
terers changes completely between successive shots. In 
this time inverval (100 ms), the scatterers diffuse a dis¬ 
tance comparable to their separation, about 1 /im. To 
remove this source of fluctuations, we arrest scatterer 
motion by gelling the liquid. A gelatinous gain medium 
allows for the free diffusion of water and dye, but blocks 
the translation of larger particles. We use a gel which is 
92% liquid and 8% polyacrylamide with a Young’s mod¬ 
ulus of about IKPa |27|. Under these conditions, the 
scattering strength and the amount of gain material in 
the gelatinous sample are nearly identical to those of the 
fluid sample, but now the scatterers are fixed with resid¬ 
ual fluctuations of about 5 nm. 

The fluctuations of emission from the gel sample are 
distinct from the liquid sample, as shown in the second 
row of Fig. The discrete lasing peaks do not change 
location from shot-to-shot. Fig. [^e), and do not aver¬ 
age into a smooth spectrum upon averaging across shots. 
Fig. [^f). In each shot, the emission is randomly dis¬ 
tributed among discrete peaks, and there is no correla¬ 
tion in single-peak height fluctuations from shot-to-shot. 
The single-wavelength intensity distributions do not show 
the slow decay characteristic of Levy statistics that was 
observed for the liquid sample. Fig. [^g). While the sta¬ 
tistical distributions of emission intensities at individual 
wavelengths in the liquid sample are very similar, the me¬ 
dian intensities vary widely across wavelength in the gel 
sample. When the emission data are pooled across wave¬ 
lengths, Fig. Bh), a power-law distribution re-emerges 
with exponent u ^ —1. However, when the intensities 
at each wavelength are normalized to the average of that 
wavelength (/(A))^ and then pooled across wavelengths, 
the distribution becomes similar to that of a single wave¬ 
length and decays much more rapidly than the power-law. 

Why do the liquid and gel gain media show different 
emission statistics? The gain materials and pump con¬ 
ditions are the same in both cases. The two samples 
have identical scatterers embedded in media with nearly 
identical optical properties. Therefore, we do not expect 
the statistical differences between these samples to re¬ 
flect differences in lasing action. Instead, the dominant 
difference between these two samples lies in the dynam¬ 
ics of the scatterers. These dynamical differences have 
no significant impact on lasing action within a single 
shot because samples have essentially stationary scatter¬ 
ers over the duration of a single pump pulse: over this 
time interval (30 ps), the particles diffuse about O.I A. 
Therefore, the remaining difference lies in the ensemble 
of scatterer configurations sampled by the series of pump 
pulses. While the particle configurations are uncorrelated 
from shot-to-shot in the liquid gain medium, they are es¬ 
sentially identical across all shots in the gelatinous gain 
medium. This explains our observation that the location 
of the peaks in the emission spectra of the gel are consis¬ 
tent from shot-to-shot. To confirm that the differences in 
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FIG. 1: Statistics of Random Lasers with Different Scatterer Motion. Data from a liquid sample, stationary gel sample and 
100 jam/s translating gel sample are shown in the first row (a-d), second row (e-g) and third row (i-1), respectively. The first 
column (a,e,i) shows emission spectra for 200 consecutive shots. Second column (b,f,j) shows the mean energy over 500 shots 
at each wavelength (thick black line), and the intensities collected in 3 individual shots (red, green, and blue lines). The third 
column (c,g,k) shows distributions of laser emission intensity for individual wavelengths between 608 and 614 nm. The colored 
distributions correspond to wavelengths highlighted by the same colored dots in the second column. The fourth column (d,h,l) 
shows the intensity distribution when data are pooled across wavelenths. The red curve shows the distribution of the intensities 
pooled across shots and wavelengths, normalized by the mean over all shots and wavelengths {I)^ x • black curve shows 
the distribution of intensities that are first normalized by the average intensity at each wavelength (/(A))^ and then pooled 
across wavelengths. 


the statistics of the liquid and gelatinous media are due 
to particle motion, we acquired random lasing spectra 
from a gel sample while translating the pump spot on the 
sample at a speed of 100 jam/s. At this speed of trans¬ 
lation, there is little overlap in the pump volume from 
shot-to-shot, so the sampled scatterer configurations are 
independent. As shown in the third row of Fig. the 
emission spectrum of translated gel recapitulates all of 
the features of the liquid sample. 

Having experimentally removed the effects of scatterer 
motion, we quantify the shot-to-shot fluctuations of indi- 
viudal lasing mode intensities. Locating peaks in the 
mean emission spectrum of the stationary gel sample 
(Fig. [^second row), we identify 25 distinct lasing modes. 
The survival functions, 5'(/, A), for each mode are shown 
as colored lines in Fig. [^. When we normalize the in¬ 
tensity of each mode by its mean and plot the resulting 
survival function, we find that all but the brightest modes 
have an exponential distribution of intensities, as shown 
in Fig. Therefore, at the level of a single lasing 


mode, we find no signature of Levy statistics. However, 
the mean intensity varies widely from mode to mode, so 
that if we pool the data across all modes before com¬ 
puting the survival function, the statistical distribution 
of emission intensity displays a slower-than-exponential 
decay, as shown by the thick black line in Fig. Hf- 
We also find that the fiuctutations of emission inten¬ 
sity are coupled across lasing modes. We quantify cross¬ 
correlations in the lasing mode intensities using the co- 
variance: 

C 2 (Ai, A 2 ) = (/(Ai)/(A 2 ))t - {I{Xi))t{I{X2))t. (1) 

The lasing modes can be significantly correlated, even 
when they are well-separated across wavelengths, as seen 
in Fig. H- Since the intensity distribution of each mode 
is nearly exponential, C 2 (A,A) ^ (/(A))|. Thus, the 
brightest modes have the highest variances, which are 
shown along the diagonal in Fig. H- These bright modes 
are anti-correlated with other modes, indicated by the 
off-diagonal elements. To scale out the trivial dependence 
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FIG. 2: Statistics of random lasing modes in a gel sample with immobile scatterers (first row, a-d) and a simple mean-field 
simulation (second row, e-h). First column (a,e): colored curves show the survival function of intensities for all 25 lasing modes 
on a log-log plot. The thick black curves show the survival functions of intensity pooled across all of the modes. Second column 
(b,f) shows survival functions of intensity normalized to the average of each mode on a semi-log plot. Third column (c,g) 
shows the correlations (covariance) across lasing modes. Red (blue) indicates positive (negative) correlations. The diagonal is 
the variance of each peak, and the off-diagonal shows correlations across peaks. Fourth column (d,g): scatter plot of Pearson 
correlations between lasing modes (excluding self-correlations). The data points with red circles indicate correlations with the 
brightest lasing mode. 


of the cross-correlations on the variance of each mode, we 
calculate the Pearson correlation between each pair of 
modes, p(Ai,A 2 ) = C 2 (Ai, A 2 )/yc 2 (Ai, Ai)c 2 (A 2 , A 2 ). As 
shown in Fig. |^, the strength of the Pearson correlation 
between modes scales with the geometric mean of their 
intensities. The strongest cross-correlations are those in¬ 
volving the brightest mode, indicated by red circles in 
Fig. iji. 

The brightest lasing modes can saturate the optical 
gain and suppress lasing in other modes, leading to the 
anti-correlations of intensities. Moreover, the intensity 
distribution of the brightest modes decay faster than ex¬ 
ponentially, as a result of gain depletion. To construct 
an empirical model to capture these features, we take 
a mean-field approach. For each pump pulse, we as¬ 
sign an random initial intensity to the n-th mode, e.g.^ 
In^. We assume that the statistical distribution of the 
initial intensity of each mode is exponential, p{In^) = 
(In exp[—/n V(^n ^)]- Since the pump pulse energy is 
identical for each shot, the total emission intensity should 
be constant. We determine the observed intensities by 
renormalizing the initial intensities of individual modes 
by the sum of the initial intensities across the modes: 
In = This normalization enforces global 

energy conservation, thus taking into account the mode 
competition for a limited gain. To compare this simple 
model to our data, we assign the experimental value of 
the observed mean intensity for each lasing mode {In)t 
to {In^). The results for this simple model with 500 
simulated shots are shown in the second row of Fig. 


The model captures the essential features of the one-point 
statistics: most of the modes decay exponentially, while 
the brightest ones decay more quickly due to gain satura¬ 
tion. Importantly, the pooled-intensity statistics across 
the modes show a slower-than-exponential decay. The 
simulation also captures the correlated intensity fluctu¬ 
ations between lasing modes, with the strength of the 
anti-correlation increasing with the geometric mean of 
the intensity of the two modes. 

This simple model faithfully reproduces the observed 
one- and two-point statistics of random lasers with fixed 
scatterers. Since a mean-field enforcement of energy 
conservation effectively simulates gain competition and 
saturation, we expect our experiment has a strong spa¬ 
tial overlap of the lasing modes. However, to specify 
the mechanism of mode-competition, the self- and cross¬ 
saturation coefficients of lasing modes will need to be 
calculated by the first-principle scattering and lasing the¬ 
ory [28l[29], taking into account their spatial and spectral 
overlaps m- A rigorous physical model remains to be 
developed to address the exponential distribution of the 
inherent intensities and the wide variation of mean in¬ 
tensities across modes. 

More generally, this work empirically separates con¬ 
tributions to random laser fluctuations from scatterer 
motion and inherent optical fluctuations such as spon¬ 
taneous emission. For the parameters probed in this ex¬ 
periment, we find that the inherent intensity distribution 
of each random laser mode is exponential, and power-law 
distributions emerge only when scatterers move or data 
are pooled across modes with widely varying average in- 
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tensities. These results provide physical insight into ran¬ 
dom laser statistics, and will stimulate further studies. 
Our results also indicate that the statistics of random 
lasing spectra reflect the microscopic motion of underly¬ 
ing scatterers. This suggests that random laser statistics 
may be able to quantify the motion of microscopic par¬ 
ticles in dense suspensions where conventional dynamic 
light scattering methods are limited. 

We acknowledge helpful discussions with Heeso Noh, 
Brandon Redding, Yindong Chong and A. Douglas Stone 
and financial support from the donors of the ACS 
Petroleum Research Fund (50872-ND6), the National 
Science Foundation (CBET-1236086) and the Office of 
Naval Research (MURI SPOOOl 135605). 


* eric.dufresne@yale.edu 

[1] A. Lagendijk, B. van Tiggelen, and D. S. Wiersma, Phys. 
Today 62 , 24 (2009). 

[2] H. Cao, Waves in random media 13, R1 (2003). 

[3] D. S. Wiersma, Nature physics 4, 359 (2008). 

[4] M. Noginov, Solid-state random lasers, vol. 105 (Springer 
Science & Business Media, 2005). 

[5] M. Noginov, S. Egarievwe, N. Noginova, H. Caulfield, 
and J. Wang, Optical Materials 12, 127 (1999). 

[6] B. Redding, M. A. Choma, and H. Cao, in SPIE 
NanoSeienee-h Engineering (International Society for 
Optics and Photonics, 2011), pp. 80950P-80950P. 

[71 H. Cao, J. Xu, E. Seelig, and R. Chang, Applied Physics 
Letters 76 , 2997 (2000). 

[8] B. Redding, M. A. Choma, and H. Cao, Nature photonics 
6 , 355 (2012). 

[9] D. Wiersma, Nature 406, 132 (2000). 

[10] N. M. Lawandy, Photonics Spectra 28 , 119 (1994). 


[11] R. C. Poison and Z. V. Vardeny, Applied physics letters 
85 , 1289 (2004). 

[12] R. Poison and Z. Vardeny, Journal of Optics 12, 024010 
( 2010 ). 

[13] H. Cao, Journal of Physics A: Mathematical and General 
38 , 10497 (2005). 

[14] K. L. van der Molen, A. P. Mosk, and A. Lagendijk, 
Physical Review A 74 , 053808 (2006). 

[15] S. Mujumdar, V. Tiirck, R. Torre, and D. S. Wiersma, 
Physical Review A 76 , 033807 (2007). 

[16] X. Wu and H. Cao, Optics letters 32, 3089 (2007). 

[17] S. Lepri, S. Cavalieri, G.-L. Oppo, and D. S. Wiersma, 
Physical Review A 75 , 063820 (2007). 

[18] X. Wu and H. Cao, Physical Review A 77, 013832 (2008). 

[19] O. Zaitsev, L. Deych, and V. Shuvayev, Physical review 
letters 102, 043906 (2009). 

[20] R. Uppu and S. Mujumdar, Optics letters 35, 2831 

( 2010 ). 

[21] R. Uppu, A. K. Tiwari, and S. Mujumdar, Optics letters 
37 , 662 (2012). 

[22] G. Zhu, L. Gu, and M. Noginov, Physical Review A 85, 
043801 (2012). 

[23] R. Uppu and S. Mujumdar, Physical Review A 87 , 
013822 (2013). 

[24] S. Lepri, Physical review letters 110, 230603 (2013). 

[25] E. Ignesti, F. Tommasi, L. Fini, S. Lepri, V. Radhalak- 
shmi, D. Wiersma, and S. Cavalieri, Physical Review A 
88, 033820 (2013). 

[26] R. Uppu and S. Mujumdar, Phys. Rev. Lett. 114 , 183903 
(2015). 

[27] Y.-l. Wang and R. J. Pelham, Methods in Enzymology 
298 , 489 (1998). 

[28] H. E. Tiireci, L. Ge, S. Rotter, and A. D. Stone, Science 
320 , 643 (2008). 

[29] H. E. Tiireci, A. D. Stone, L. Ge, S. Rotter, and R. J. 
Tandy, Nonlinearity 22, Cl (2009). 



